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Abstract. We investigate the three-dimensional formulation of a recently proposed 
operational arrival-time model. It is shown that within typical conditions for optical 
transitions the results of the simple one-dimensional version are generally valid. 
Differences that may occur are consequences of Dopplcr and momentum-transfer 
effects. Ways to minimize these are discussed. 
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1. Introduction 

In quantum mechanics, as opposed to classical physics, the meaning of arrival time of a 
particle at a given location is not evident when the finite extent of the wave function and 
its spreading become relevant, such as for cooled atoms dropping out of a trap. Moreover, 
in the quantum case one will expect an arrival-time distribution, and there are different 
theoretical proposals for it, see e.g. [1112113111001111313110111111121, the review 
in ^3], and the book ^3]. These arrival-time distributions have been controversially 
discussed in the literature since they are derived through purely theoretical arguments 
rather than operationally, i.e. without specifying a measurement procedure. But from 
the experimental point of view this is of course very important. 

How then to measure the arrival time of a single quantum mechanical particle at a 
given location when the extent of its wave function becomes important? Here one has 
to distinguish between first arrival and repeated arrivals as in the case of a pendulum. 
For the former, the measurement procedure has to ensure that the arrival is indeed the 
first, not the second or third, i.e. that the particle did not already arrive before. This 
implies a sort of continuous or rapidly repeated observation of the same single particle 
in question. Doing the experiment again and again with many identically prepared 
particles will then give an arrival-time distribution. 

To capture these requirements in a model it was suggested in [12] to consider a 
particle ( "atom" ) with two internal levels, which are connected by an optical transition, 
and to illuminate the arrival region by a laser. When the wave function of the particle 
enters the region the laser will populate the upper level and one can watch for the first 
spontaneous photon. The detection time of this photon can then be taken as the first 
arrival time in this region, at least approximately. 

The model contains of course simplifying idealizations, which are indeed quite 
essential for extracting useful information. In general a measurement model should be 
flexible enough to capture the essence of the experiment and still provide, hopefully, 
close links with the schema of fundamental theory. That the fluorescence-based 
model is of such kind has been already demonstrated: the limits or operations by 
which theoretical time-of-arrival distributions, such as the flux or current density and 
Kijowski's distribution, can be obtained by measurement have been determined [TKllKij. 
For other applications or extensions of the model see [T71 [13 113 123 121] • It should not be 
surprising that the conditions required for these "perfect measurements" may be rather 
extreme, or complicated to realize. The important point is that at least one knows what 
should be done, and one can reasonably understand and try to minimize the deviations 
from the ideal results. 

A very basic simplifying assumption of all applications of the model so far has been 
the restriction of the atomic motion to one dimension (ID), associated with the x axis 
hereafter, perpendicular to the laser (traveling wave) direction y. This approximation 
could in principle be badly wrong because of transversal velocity components in the 
incident atomic wave function or because of momentum transfer from the laser. 
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Figure 1. Schematic picture of the atom wave packet moving towards the laser- 
illuminated region. 



The aim of this paper is to generahze the fluorescence model to an arbitrary 
three dimensional wave packet impinging on a plane (x = 0) which separates a laser- 
illuminated region [x > 0) from a non-illuminated region (a; < 0). It will be shown that 
possible deviations from the ID model can be attributed to an additional detuning term 
which is of kinetic origin. This kinetic detuning consists of a Doppler term and a term 
which is due to a momentum transfer from the laser field to the atom. It will be seen that 
deviations of the time-of-arrival distribution from the ID model can be ascribed to two 
effects. First, there is increased or decreased reflection for blue or red kinetic detuning, 
respectively and, second, there is a less efficient driving of the atomic transition for 
both signs of the kinetic detuning which leads to a delay in the distribution. Numerical 
examples will be shown for two typical cases. For oblique incidence the kinetic effects 
can be compensated by appropriate laser detuning, whereas for a wave packet prepared 
with a very small width in laser direction this is not possible. The numerical results 
also show that the ID model is valid over a wide range of parameters. 

It will also be shown that the first photon distribution tends in some limit to the x 
component of the total fiux through the plane x = 0, which is the natural generalization 
of a main result in [TK] . 

2. The 3D model 

As in we consider a two-level atom which moves from the left towards the region 
X > which is illuminated by a traveling-wave laser directed parallel to the y axis (see 
figure P). The laser frequency may be on resonance or be detuned with respect to the 
atomic transition. We use the dipole and rotating wave approximations and employ the 
quantum jump approach |22l EHl Elj • In this approach the time development between 
two photon detections is given by a nonhermitean "conditional" Hamiltonian so that 
the norm of the wave function decreases; the norm squared at time t just gives the 
probability for no photon detection until t. In a laser-adapted interaction picture (i.e. 
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with Hq — hu!i^\2){2\) the conditional Hamiltonian becomes 

^ + + v(±) (1) 

' 2m 2m 2m ^ ' ^ ^ 

with 

V{±) = ln9{x) [e''=i^^|2)(l|+e-''=i^*|l)(2|l -;i(AL + i7/2)|2)(2|, (2) 

where 7 is the decay (Einstein) constant of the excited state |2), 1] is the Rabi frequency, 
Al the detuning (laser angular frequency minus atomic-transition angular frequency), 
/cl the laser wavenumber, and 9 the Heaviside function. The hats on positions and 
momenta (x, x, y, ...,Px, ■■■) denote operators to be distinguished from the corresponding 
ordinary vectors or c-numbers. Boldface capital Greek letters will designate atomic two- 
component wave functions, with the components for excited and ground state denoted 
by superscripts. 

The operational time-of- arrival distribution n(t) for an atomic ensemble represented 
by a given wave function is defined as the temporal distribution of the first photon 
detected for each atom, and this is just the rate of the decrease of probability of detecting 
no photon. Hence one has the general result 

n(t) = -d||*(t)||Vdt . (3) 

By means of the Schrodinger equation and the above form of He one easily finds for the 
two-level system 

-00 

To calculate H it is useful and physically illuminating to first obtain the eigenstates of 
the conditional Hamiltonian, 

i7,#k = EA- (5) 

Since K(x) does not depend on z, the z component of momentum is conserved so that 
there is free motion in 2;-direction and we only have to consider the two-dimensional 
motion in the x-y-plane. Therefore, from now on the position and momentum 
(wavenumber) vectors will always be two-dimensional. 

Since the incident atoms are initially in the ground state and far from the laser 
region and are moving in the positive x-direction, we need the stationary scattering 
eigenstates which correspond to ground state plane waves coming in from the left. 
In matrix form, with |1) = and |2) = the (two-dimensional) conditional 
Hamiltonian can be written as 

Pi Pi • / o\ h^^^^f ne-'''^y \ 

^ 2m 2m ^ ^' ' \ 1 2 ^ ' \ Oe'^^^y / ^ ^ 



For X < we use the ansatz 




(7) 
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with (27r) ^ included for delta normalization in two-dimensional k space and with 

= ^ = ^ q^ = k^ + -(i-, + 2A0. (8) 

At this point the individual components of k' and q are still unknown. They will be 
determined later from the matching conditions at x = 0. For the region a; > we use, 
in analogy to the ID model [T^ I18j . a plane wave ansatz of the form 



^k(x)=( ^A,^^,,,, ]^^^-^\E^y^^- (9) 



with as yet unknown A± and k^. Inserting this into the eigenvalue equation (0) the 
exponential e''^''''^ drops out and we get the matrix equation 



— k^^ -Vte~'^^^y 

2m 2 

pe-^'^y |ik±2 - ^(A + i7/2) 



|E±) = Ek|S±) , (10) 



where 



A ^ - ^{2kp^ + kl) (11) 

can be regarded as an effective detuning. We shall later analyze in detail the physical 
meaning of the second term, which is clearly of kinematic nature. Equation (|1U|) is 
fulfilled if X± and k^ satisfy 



A± = (i7 + 2A) ± -V(7 - 2iA)^ - 4^^ (12) 

and 

h 

Note that X± still contains the unknown components fc^ and k^. The eigenstate for 
X > is now a superposition of the form 

$ij(x) = ^ [C+|E+)e''^^-" + C_|^_)e''^"-"1 . (14) 
The coefficients C± will be determined from the matching conditions. 

2.1. Matching Conditions 

At X = the usual matching conditions, i.e., 

<(x = 0,y) = $L''(a^ = 0,2/) 
are imposed. This yields the equations 

e'''yy + Rie'^> = C+e'''yy + C^e^'^y^ (16a) 
R^e^'^yy = ^C'+e'('=L+4)j/ + ^C'_e^('^L+fc„-)s/ (Ig^) 

ke''^^^ + k'Rie'''yy = k+C+e^'^^^^' + k-C^e^'^y^ (16c) 
qR^e^'^yy = k+^C+e^('=L+fc+)y + i^-'^c^e'^'^^+'^y^y (16c/) 



(15) 
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As (I16aj) has to hold for all y one obtains 

= = ky = ky , (17) 
and with this one finds from (jHJ for k'^ of the refiected component 
k^ = —kx- 

From (fTTjl it follows that A is independent of k^, and consequences of this will be 
discussed in more detail in the next subsection. For the x components of the wave 
vector inside the laser region one then obtains 

m , , , im 



fcf = + ^ (i7 + 2A) T §f\/(7 - 2iA)^ - m\ 

and fc^ are given by the roots of this with positive imaginary part so that the wave 
decays in the laser region. Both depend on ky through the effective detuning A. From 
()16?)|) one gets by a similar argument 

Qy — ky -\- k]^, 

where hki^ is the momentum transfer from the laser field to the excited state of the 
atom. This and ^ in turn lead to 

ql = kl + j{i^ + 2A) 

for the refiected part of the excited state, and is obtained by again taking the root 
with positive imaginary part so that the excited wave decays at long distances from the 
laser. With these simplifications the matching conditions (16) now take the form 

1 + R^ = C+ + C^ (18a) 

kxil-Ri) = k+C+ + k-C. (18c) 
2A 2A 

which look formally just same as in the ID case, only with k^ and q replaced by k^ 
and qx and with A± defined as in (fT^. Accordingly the coefficients Ri, R2, C± have the 
same form as in the ID model, 

C+ = -2kx{qx + k-)\-/D (19a) 

C^ = 2kx{qx + kt)\+/D (19&) 

i?2 = kx{k- - k^)n/D (19c) 

Ri = [A+(g. + kt){kx - k-) - A4g. + k-){kx - k+)]/D (19c?) 

and 

D = A+(A;, + k-){qx + kl) - X_{kx + kt){qx + k'). (19e) 
Finally, taking into account the matching conditions, ((Tj) and (|14p take the form 




gi/c^z + R^e 



R2e 



-igxx i{ky+kL)y 



(20) 
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for X < and 




1 




+ a 



1 

n ^ 



(21) 



for X > 0. 



2.2. Effective detuning: Doppler shift and momentum transfer 



Using (HZj), the effective detuning (fTT|) can be written as 

A = Al - Ak, 
with the kinetic detuning 



(22) 



SO that ^Ak is the y momentum component contribution to the kinetic energy difference 
between ground state and excited state, as seen from the terms ky and ky + /cl in the 
plane wave exponents in (j^UI) and (j?T|) . We may also express Ak in terms of the speed 
of light in vacuum c, the laser angular frequency cjl = /clc, the y component of the 
incident velocity Vy = hky/m, and of a velocity vn = hki^/m, namely 



The form of ()23|1 can also be derived by classically modeling the resonance condition 
with simple energy and momentum conservation arguments. To see this we rewrite it 

as 



with ujd = uoi^Vy/c and cjr = hk1/2m = uji,vyi/2c. The first term is the ordinary Doppler 
shift due to the incident velocity component Vy and the second can be viewed as a 
frequency shift caused by the change of the atomic kinetic energy due to the absorption 
of a laser photon. Recall that the momentum transfer from the laser is only directed 
along the y direction. It is interesting to note that the kinetic detuning in ()24|) has in 
fact already the form of a Doppler shift formula with the role of the velocity played 
by the average of the velocity components for ground state and excited state along the 
direction of the laser. Indeed, for fixed non-zero values of ky and /cl in a stationary 
wave the two effects cannot be disentangled and act simultaneously in the form of a 
single kinetic detuning. It is, however, possible to distinguish their consequences for 
wave packets, as discussed below, since is fixed by the laser wavelength, whereas ky 
varies for each incident atomic momentum component, or by sweeping over the incidence 
angle. 

Note also the signs in the effective total detuning ()22p : a large, positive, kinetic 
detuning (energy transfer) implies a red shift whereas a large laser, positive detuning 
amounts to a blue shift. They may compensate each other, as we shall illustrate 
afterwards, while and cjr may also cancel each other for a negative velocity Vy 
opposed to the laser beam. 




(23) 



Ak = cjd + i^R 



(24) 
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3. The First-Photon Distribution 



Until detection of the first photon, the conditional time development of a wave packet 
^'(x, t) corresponding to a particle in the internal ground state coming in from the left 
can now be written as in the ID case in terms of the plane-wave solutions, 

dA;,^(k)$k(x)e-'"'^'*/2'", 



*(x,t) 



(25) 



where ip(k) is the momentum amplitude the wave packet would have as a freely moving 
particle at t = 0. Inserting this into ^ one obtains the first-photon distribution. This 
gives a six- dimensional integral of which the x, y, and k'y integrations can be carried out 
analytically, yielding 



n(t) 



/y roo poo 

Zir Jo Jo 



dky ip[kx, ky)ip[k^, ky 



ii?2^2 



+ i- 



AC+C X+X' AC-C'X^X' 



ACZC'_X^X'_ 



■ih{k 



(26) 



X 

'2_i.2 



ki)t/2m 



^i\k-' - kt) n^ikr - k-) n^{k-' - k-,. 

This form of the first-photon distribution will be now used to numerically look for 
deviations from the ID results in different parameter regimes. 

The relevance of possible deviations can be estimated by comparing the two- 
dimensional (2D) eigenstates with those of the ID model which are given by IT8] 

-ikx \ 



$r(x) 



where 











1 


L V f7 


V2^ 



+ c_ 



2A 



(ID) 



Ak-X 



for X < 



for X > 



(27) 



/p + ^(AL + i7/2) 



A. 



-- (i7 + 2Al) ±^^J{l- 2iAL)2 - 4fi2 



and roots are taken with positive imaginary parts. One can go from 2D to ID formulas 
{Qx X^2^^ ~^ ^±^\ ~^ k±) simply by dropping the kinetic detuning Ak, which 
is thus identified as the physical reason for possible deviations from the ID results. 

For optical transitions, typical values of Q, 7 or Al are of order 10'' s^^, while the 
shift ujn in ()23p for Cs is of order 10^ s~^. For perpendicular incidence, i.e. (ky) = 0, 
and for velocities Vy corresponding to a Ay = 1/im wave packet, the Doppler shift o^d is 
even smaller, namely of order 10^ s~^. So the ID model results may be expected to hold 
quite generally in a broad (ordinary) range of parameters. The number of parameters 
is rather large so a full systematic analysis of all possible cases in the vast parameter 
space is out of the question. Nevertheless, it is worth examining the convergence from 
2D to ID results in some typical cases. Since the laser and kinetic detunings, Al and 
Ak, are always combined together in X±, fc^ and q, one possible limit for pure ID 
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Figure 2. Arrival time density for a caesium atom with the parameters f7 = 
1.67 • 10^ s~^, 7 = 3.3 • 10^ s~^, prepared at i = as a minimum-uncertainty product 
Gaussian wave packet with = 0.24 • 10^^ m, {k^) ~ 0, and with initial average 
position (xq) = —1.32 • 10^^ m, and velocity (vq) = 9cm/s. The solid line distribution 
coincides with the ID distribution. 

behaviour is |Al| ^ |^k|- Moreover, from ()22|) it is also possible to find ID results 
if the two contributions to the effective detuning cancel each other, i.e. if Al = Ak, 
as we shall discuss below. Ak may also be zero because of a cancellation between 
Doppler and momentum-transfer shifts in oblique incidence. One may similarly reduce 
the equations to ID if f2 or 7 are large with respect to all other frequency parameters. 
For strong driving, 7, Al, more detailed sufficient conditions would be ky <C k^/ki^, 
ky,ki^ -C ^Q/ki^, and kf^ <^ k^. This means that a two-dimensional wave packet, with 
given mean momenta (k^), (ky) and momentum widths A^^, A^^^, can be described by 
the ID model if it has only small ky components and if at the same time /cl is small 
compared to the k^ component and Ak ^ In particular the momentum spread of 
the wave-packet in y direction should be small compared to the momentum spread in x 
direction. This, on the other hand, corresponds to a large width Ay in real space which 
is very satisfactory since for Aj^ 00 the two-dimensional wave packet tends to a one 
dimensional one. 

This behaviour can be seen in figure El for arrival time distributions of minimum- 
uncertainty product Gaussian wave packets with different widths A^ and fixed width 
Ax at the preparation time t = (see the caption for details), where parameters of Cs 
have been used. For large widths Ay the time of arrival distribution is identical to the 
one which is obtained from a one-dimensional Gaussian wave packet of width A^; in the 
ID model. For very small widths Ay first a slight enhancement in the height of the 
distribution and then a delay in the form of an enhanced tail of the distribution can 
be seen. The physical reasons for these two effects are a dependence of the reflection 
coefficient on the velocity in y direction and the Doppler effect, respectively, as will 




— 
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Figure 3. Arrival time density for Gaussian wave packets with different positive 
initial average velocities (vy). Atomic parameters as in figure |21 The dotted graph 
shows a distribution with Al = Ak for compensation of the kinetic effects. 



be explained in the next paragraph. It is seen in figure |21 that a significant deviation 
from the ID distribution occurs for the Cs parameters used here if the wave packet is 
prepared with a width at least three orders of magnitude smaller in y direction than in x 
direction. A deviation from the ID time of arrival can also be seen if a sufficiently large 
momentum (ky) is chosen, which corresponds to oblique incidence of the atoms onto the 
laser region. Figure 01 shows numerical results for Cs wave packets with different positive 
average initial momenta. Again there are two effects: first an increase in height of the 
distribution and then for even larger momenta a delay. For negative initial average 
momenta (see figure E} on the other hand one first sees a decrease in height of the 
distributions and then also a delay. 

The change in height of the distributions is due to the ky dependence of the refiection 
coefficients, as will be explained now. Compared to the eigenstate with ky = the 
eigenfunctions with positive ky have a smaller refiected part whereas the eigenfunctions 
with negative ky have a larger refiected part, as shown in figure El This behaviour 
can be understood via the dipole force and the effective detuning of (j^^ as follows. 
Positive values of ky lead to an effective red detuning while negative values lead to blue 
detuning. Therefore states with positive ky are pushed into the laser region resulting 
in decreased reflection and states with negative ky are pulled out of the laser region 
resulting in increased reflection. The delay is due to the diminished efficiency of the 
laser because of the Doppler shift which also explains the decrease in refiection for very 
high negative momenta seen in figure El For the parameter values of Caesium used 
here one again has to insert rather extreme values of (ky), i.e. at least three orders of 
magnitude larger than (k^), to see deviations from the ID distribution. The effects of 
nonzero initial average momentum ky can be compensated by choosing a laser detuning 
Al = Ak which leads to zero average effective detuning. Since quantum effects in the 
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Figure 4. Arrival time distribution for wave packets with different negative initial 
average velocities (vy). Atomic parameters as in figure |21 The dotted graph again 
shows a distribution with Al — Ak for compensation of the kinetic effects. 



time-of-arrival can only be observed for sufficiently small incident velocities and since 
the 3D distribution corresponds to the ID distribution with as incident velocity the 
oblique incidence could be useful in order to obtain broader distributions. 

The effects of very small widths Aj^, seen in figure |2l can also be understood by 
means of the two effects explained above. The delay is again due to the Doppler effect. 
Indeed, since small uncertainty in position means large uncertainty in momentum, a 
wave packet prepared with small will broaden very rapidly so that the tails of the 
wave packet move away from the center. Due to the Doppler effect these tails see a 
detuned laser and the wave packet as a whole is excited less efficiently. Since there is 
a shift with positive and negative sign for different components, it is not possible to 
compensate this effect by laser detuning. The increase in height of the distribution is 
a little bit subtler. Starting from ky = the decrease of in figure El when 

going to positive values of ky is slightly steeper than the increase in the other direction. 
Therefore there is less reflection for increasing Aky, although both positive and negative 
values of ky occur with the same weight in the wave packet. 

Deviations from the ID distribution also occur for a large shift wr, i.e. if k^ is large 
compared to the other relevant parameters. However, this requires an incident velocity 
smaller than the recoil velocity vn and, for a laser transition in the optical range, a 
metastable transition with a lifetime of the order of one tenth of a second. 

4. Deconvolution 

The operational model to quantum arrival times presented in and generalized here 
to three dimensions shows that there may arise two potential problems when comparing 
operational with ideal distributions, namely reflection due to the laser field and detection 
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delay. Reflection in the ground state means that no photon is emitted and hence the 
arrival of the atom is not detected. Delay is due to the fact that excitation and de- 
excitation of the atom take a finite time. If one tries to get rid of the delay by increasing 
both the decay rate 7 and the Rabi frequency Q this will increase reflection in the ground 
state and thus lead to a detection decrease. To circumvent this one may try to use a 
gentle (i.e., not reflecting) weak excitation and then subtract the long delays in a suitable 
manner. In [inj this was done by a deconvolution with the flrst-photon distribution of 
an atom at rest as follows. It was assumed that the "experimental" time-of-arrival 
distribution 11 (t) is given by a convolution of a hypothetical ideal distribution Hid with 
the (known) flrst-photon distribution W{t) of a two-level atom at rest, 

n = Hid * W. (28) 

In the limit of large 7 it was then seen that the resulting ideal distribution tends to the 
quantum mechanical flux 

Hid -^J^ = ^^ {mt)i^'io, t) - WiMm t) } . 

A similar calculation can be carried out here and it yields the x component of the total 
flux through the plane x = 0, 



Hid ^ Jx(0,t) 

with 



dyJ^{x,y,t) 



x=0 



-dip dip 



Jx{.x, y, t) = ^ ^ip{x, y, t)-^ix^ i) - ■^(^' ^' y^t)j. 

This is the result one would expect as the natural generalization when carrying over the 
ID case to three dimensions. 



Quantum optical time- of- arrival model in three dimensions 



13 



5. Conclusions 

We have generalized the quantum optical time-of-arrival model of ^H] to three 
dimensional space in order to describe a more realistic setup and verify the 
approximations made there. When considering arrivals at the plane x = we have 
seen that only x and y directions have to be considered, which leads to an effective 
2D model. Deviations from the ID model can be described via a kinetic detuning Ak 
which consists of two terms for two distinct effects. One is a detuning of the atomic 
transition from the laser frequency due to the Doppler effect for velocity components Vy 
in the laser direction. The other is the gain of momentum of the atom in the positive-?/ 
(laser) direction by absorption of a laser photon. Again via the Doppler effect, this 
leads to a shift in frequency of the internal transition of the atom in the lab frame. 
A nonzero kinetic detuning in turn leads to two different effects affecting the time-of- 
arrival distribution. Blue and red detuning results in increased and decreased reflection, 
respectively, and therefore in a smaller or larger height of the distribution. At the same 
time both detunings lead to a less efficient driving of the atomic transition by the laser 
and thus to a delay in the distribution. 

We have described typical situations for which deviations from the ID model 
occur. One is the preparation of a wave packet with very small width Aj^. Then the 
corresponding large momentum width leads to large transversal velocity components. 
Another possibility is oblique incidence with large (positive or negative) mean 
transversal momentum (ky). In the latter case the deviations can be compensated by an 
appropriate detuning of the laser whereas in the first case this is not possible. Inserting 
realistic atomic parameters it has been shown that over a wide range of parameters the 
ID model is generally applicable. 

A main theoretical result of reference J3] was that the first-photon distribution 
tends in some limit to the quantum mechanical probability current, opening the way 
towards a measurement of this quantity. A similar calculation in the generalized case 
gives the x component of the total flux through the plane x = 0, which is the natural 
generalization of the ID result. 
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